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Abstract - In this research work, we propose a system of nonlinear ordinary differential equation used to model the 
interaction between Langerhans cells and HIV infection. The model consists of five compartments, namely, susceptible 
Langerhans cells, infected Langerhans cells, susceptible T-cells, infected T-cells and free HIV particles. The biology of 
interactions of Langerhans cells with HIV and mathematical preliminaries which plays a crucial role in our entire research work 
was described. By presenting a theoretical framework related to the infection mechanism, a biologically meaning full 
assumption was considered. Furthermore, the positivity of the model solution, the equilibrium point (both virus free and endemic 
equilibrium) of the model was shown and its stability was investigated. Finally, by using a numerical simulation the developed 
model was studied and the results concluded that the numerical simulation matches the analytical solution as expected. 
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1. INTRODUCTION 


The term Langerhans cell (LC) was studied by a German medical student Paul Langerhans. He was 
21 years old when he described in 1868 a new epidermal cell type which now bears his name. The cells 
have long dendrites, so Langerhans mistakenly regarded them as a component of the nervous system. A 
century later, Inga Silberberg discovered that these cells instead have a role in immunity, and later work 
showed they are a type of dendritic cell, a classic antigen-presenting cell. A network of epidermal LC 
laces through the skin and mucosa, including the anal and vaginal mucosa as well as the male foreskin 
[2]. 

An electron microscopic study showed that during sensitization to allergenic compounds, LC 
migrate from the skin through the lymph vessels and come into close contact with lymphocytes in the 
draining lymph nodes [3]. These findings suggested that LC play an important role in the development 
of the immune response to skin antigens and triggered various groups to investigate the possibility that 
it may function as the most peripheral arm of the immune system. Throughout the 1970s and early 
1980s, it was established that the cell represented an essential element of the immune system playing 
a central role in the mechanisms of defense of the body in a wide range of pathological processes [4]. 

It is a well known fact that the CD4+ T-cells are targets of HIV and are also important for the 
establishment and maintenance of an adaptive immune response. CD8+ T-cells are the primary 
effectors cells in HIV infection, as they kill infected cells and produce non-lytic antiviral factors. In lymph 
node, LC activates CD4+ and CD8+ T-cells during its presentation of antigen. Originated from the bone 
marrow, LCs migrate to the peripheral epithelia (skin, mucous membrane) where they play a primordial 
role in the induction of an immune response and are especially active in stimulating naive T 
lymphocytes in the primary response through a specific cooperation with CD4+ T cells after migration 
to proximal lymph node. Also LC can repopulate from local and blood-borne precursors. Similar to other 
myeloid DC, LCs is able to bridge innate and adaptive immunity. They are able to interact directly with 
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microorganisms at the periphery to produce effectors cytokines and initiate or re-stimulate activation 
of T and B cells through antigen presentation [7, 8]. 

The immunocompetent cells which act as antigen-presenting cells, Langerhans cells, can be 
infected. Its infection by HIV is relevant to several reasons. Firstly, LCs of mucosal epithelia may be 
among the first cells to be infected following mucosal HIV expositor and, secondly, LCs may serve as 
reservoir for continued infection of CD4+ T cells, especially in lymph node where epidermal LCs migrate 
following antigenic activation [5].Additionally, many indirect and/or direct experimental data have 
shown that LCs may be a privileged target, reservoir, and vector of dissemination for the HIV for the 
inoculation sites (mucosa) to lymph nodes where the emigrated infected LCs could infect T 
lymphocytes. Apart from many plasma membrane determinants, LCs also expresses CD4+ molecules 
which make them susceptible targets and reservoirs for HIV. Once infected these cells due to their 
localization in areas at risk (skin, mucous membranes), their capacity to migrate from the epidermal 
compartment to lymph nodes, and their ability to support viral replication without major cytopathic 
effects could play a role of vector in the dissemination of virus from the site of inoculation to the lymph 
nodes and thereby contribute to the infection of T lymphocytes [6]. 

Infection of biopsies of human cervical and skin of primate foreskin tissue explants show that LC 
can be infected. Topical infection of human vaginal epithelial explants with HIV strongly suggested that 
LC and non-activated T-lymphocytes are the major cell types expressing HIV antigen in and emigrating 
from the explants and they are often associated during emigration with HIV antigen concentrated at 
their contact region. The latter suggests that LCs is transferring HIV to CD4 T cells. LC may also provided 
a mode of intracellular storage while transporting HIV to CD4 T cells in the sub mucosal lymphoid tissue 
and thence to draining lymph nodes [3]. 

Alternatively, another data suggest that LC has an anti-viral function by capturing HIV for 
degradation and thus initially impairing HIV transmission. This anti-viral function is dependent on viral 
load and LC phenotype, strongly suggesting that the role of LC in HIV transmission is also likely to be 
influenced by local conditions such as viral load, stage of the menstrual cycle, state of vaginal mucus 
and/or inflammation and co-infections. They are also of particular importance because HIV exploits it 
to enhance infection. Therefore, LC are the critical link between virus, CD4+ and CD8+ T-cells. After 
encountering antigen in the periphery, LC mature and travel to the lymph node (LN). LC maturation 
includes increasing antigen presentation on major histocompatibility complex (MHC) molecules, and 
up regulating co-stimulatory molecules.Mature LC prime CD4+ T-cells to become effector T-helper cells. 
Additionally, LC cross-present exogenous antigens on MHC to prime CD8+ T-cells differentiate into 
cytotoxic T lymphocytes (CTL). Thus it is essential for fighting intracellular pathogens like HIV. In 
contrast, LC may act as one of the primary initial targets for HIV infection. Since it is specialized by 
antigen presentation and belongs to the skin immune system, the virus can associated with it to travel 
to lymphoid tissue, where 98% of T-cells reside. During antigen presentation, it can facilitate infection 
of CD4+ T-cell. As part of the normal immune response, LCs captures virions at the site of transmissions 
in the mucosa (peripheral tissues) and migrates to the lymphoid tissue where they present to naive T 
cells and hence are responsible for large-scale infections of the CD4+ T cells. Generally, as a result of 
these roles of LC, in HIV infections, it plays a dual role of promoting immunity while also facilitating 
infections [2, 8, 9]. 
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1,1. Mathematical Preliminaries 
In this section we present some definitions and theorems required to analyze model systems in our 
research work. 


Theorem 1.1 (Existence and Uniqueness Theorem): Assume D is an open subset of R X IR”, and 
f:D — R?” is continuous. Then for each ty € Rand x, € R”, the initial value problem (IVP) 
6 Gn wet a) — cn (1) 


has a solution x. If in addition, f has continuous first order partial derivative with respect to 
X4,X2,++,Xn in R”, then equation (1) has a unique solution. 

The proof is found on [10]. 

Definition 1.1: Let f is a real valued function defined on a domain D. The function f is said to be 
bounded on D if and only if there is a positive number M such that |f(x, y)| < M for all (x,y) € D. 





Definition 1.2: A domain 0 © R" is said to be positively invariant for the system = = f(x) if for 


all x(t) € 0, then x(t) € 0 forallt > to. 
Routh-Hurwitz Criteria 
Given a polynomial, 
P(A) =A" +.aya™™1 + agA”™™ +e + An_-1 + Ay 
where the coefficients a; are real constants, for i = 1,2, ...,n, define the n Hurwitz matrices using the 
coefficients a; of the characteristic polynomial: 





a, 1 00 : O 
a 1 ay 1 0 a3 az ay 0 aa 0 
H, =[a,], Hz = be a | 4s =|43 @2 a}, and Hy, = ds GQ, G3 az “+ O 
3 a2 ae a 
as Aq, a3 : i : w O 
0 0 0 0 iz. Oy 


where a; = Oif j > 0. 

All of the roots of the polynomial P(A) are negative or have negative real part if and only if the 
determinants of all Hurwitz matrices are positive: det H; > 0 for all j = 1,2,...,n. 

When n = 2, we have P,(A) = A + a,A +. a,. Hence by Routh-Hurwitz criteria we have detH, = 
a, > Oand 


1 
det H, = det a | =a,a,>0 


Oj a, 
Similarly, whenn = 5: 


. a ,az-a a1(a,a4-a 
a; > 0 fori = 1,2,...,5,2 = > 0,a, — 21(41d4=as) SQ a4a, — ag — ——— +2 _ 
a4 a4a2-a3 a3 (A1a2—a3)—a4(a1a4—-as) 


2 
as (a1a2—a3) 


>0 (2) 


2. MODEL FORMULATION AND ANALYSIS 

Aspects of an organism's defense against viral and bacterial infections and the reaction of immune 
system to infection are the main problems in practical immunology. In addition to antiviral and 
antibacterial defense, the immune system plays a decisive role in tissue incompatibility reactions, 
antitumor immunity, autoimmune diseases, and allergies [6]. Therefore, the formulation of 
mathematical models of the interaction between Langerhans cells and the spread of HIV infection 
provides the tangible understanding on these systems and the HIV transmission with its life-span. 
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Model Formulation 


In this section, we develop a model which describes the interaction between HIV and Langerhans 
cells as a system of nonlinear ordinary differential equation. For the model we denote by (L) the 
concentration of healthy Langerhans cells (or susceptible Langerhans cells), (L;) the concentration of 
infected Langerhans cells, (T) the concentration of healthy T-cells (or susceptible T-cells), (T,) the 
concentrations of infected T-cells and (V) the concentration of free HIV. 


2.1.1. Assumptions of the Model 
The main feature of the model is that the force of infection is obtained by averaging the probability 
of exposure of Langerhans cells and T-cells to HIV. For the model, the following assumptions have been 


taken: 


(i) 


(ii) 


(iii) 


(iv) 


(v) 


(vi) 


The amount of healthy Langerhans cells can increase due to newly recruited (produced) 
individuals. In converse, the number can decrease due to the natural death rate and infection 
rate. 

The susceptible Langerhans cells become infected Langerhans cells when in contact with the 
free virus and it remove from the circulation by both natural death and virus induced death 
rate. 

The number of healthy (susceptible) T-cells increases due to its constant rate of production 
in bone marrow, but the number can decreases due to its natural death rate, infection rate 
by free virus, infection rate by infected Langerhans cells and infection rate by infected T-cells. 
The susceptible T-cells become infected T-cells when in contact with free virus and infected 
Langerhans cells; and the infected T-cells remove by its natural death rate and virus induced 
death rate. 

The effective antiviral immune response of healthy T-cells in collaborations of healthy 
Langerhans cell depends on the amount of virus present, the infected tissues and the 
chronicity of the infection. 

An infected Langerhans cells and infected T-cells release hundreds of virions that can infect 
a healthy Langerhans cells and healthy T-cells. 


The schematic diagram used for the development of model is shown in Figure (1). 
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Ty 





B,LV 








pL; oT; 
Figure 1: Model structure 
Parameter description Symbol 
The constant rate of production of susceptible Langerhans cells Ty, 
from the bone marrow 
Infection rate of susceptible Langerhans cells by free virus By 
Natural death rate of Langerhans cells L 
Death rate of L; due to virus (virus induced death rate) 6 
The constant rate of production of susceptible T-cells from the Tr 
bone marrow 
Natural death rate of T cells o 
Infection rate of T-cells by free virus B2 
Infection rate of T-cells by infected Langerhans cells B3 
Infection rate of T-cells by infected T-cells Bo 
Death rate of infected T-cells due to virus (virus induced death) p 
Removal rate of free virus a 
Number of the virus particles assumed to produced by the infected N 
Langerhans cells during its life time including any of its daughter 
cells 
Number of the virus particles assumed to produced by the infected M 


T-cells during its life time including any of its daughter cells 


Table 1: The parameter description and its symbol 
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Depending on the assumptions stated and the schematic diagram of our model structure, we can 
formulate the governing model equations: 


dL 
dé — B,VL — pL 
dL, 
de = B,VL — pL, — 6L, 
aT 
ae = 1, — oT — BVT — B3L,T — B,T,T (3) 
aT, 
ae = BVT + B3L,T + B4T,T — oT; a pT; 
dV 
ae: = N6L, + MoT, = av 
where ~ — ~ at and ~ denote the rates of change of population densities of concentration of 


susceptible Langerhans cell (L), infected Langerhans cell (L,;), susceptible T-cells (T), infected T- cells 
(T;) and free HIV (V) at a time t respectively. The descriptions of symbols of parameters are explained 
in Table 1. We also suppose that w = 62V + £3L, + 6,T, for the schematic diagram shown on Figure 1. 


Zadn Model Analysis 

In this section we study the basic properties of the solutions of system of equation (3). For this, we 
are going to show the existence and uniqueness of the solution, positivity of the solution, steady states 
and stability of the solution. 


2.2.1. Existence and Uniqueness 
Theorem 2.1: Suppose that a system of equation (3) is given with non negative initial value 





X (to) = (Lo, Li, To, Ti,» Vo). Then it has a unique solution. 


Proof: We proof Theorem 2.1 by using Theorem 1.1. 
(i) Existence: 
Let (t, L, L,,T,T;,V) C R x R°. Assume 
f, CL, L,,T,T,,V) = 1, — B,VL — wh 
fo(L, L,,T,T),V) = ByVL — pL, — 6L, 
f3(L,L,,T,T,,V) = ar — oT — B,VT — B3L,T — ByT)T 
f,(L,L,,T,T,,V) = B,VT + B3L,;T + B4T,T — oT; — pT; and 
fs(L, L,,7,T;,V) = N6L; + MpT; — aV 


Then all f; are continuous in (t, L,L,,T,T,;,V) C Rx R® fori =1.,...,5. By Theorem (1.1) for each ty € 
R and (L,L,,T,T;,V) € R° there exist a solution of the system of equation (3) with the given initial 
values. 

(ii) | Uniqueness: 

Since all f; has a continuous first order partial derivatives with respect to L,L,,T,T; and V, by 
Theorem (1.1) the system of equation (3) with the given initial value has a unique solution. 
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2.2.2. Invariant Region 

Since we are modeling interaction between Pathogens and healthy cells at cellular level, the 
variables and parameters are assumed to be positive for allt => 0. The system of equation (3) will 
therefore be analyzed in a suitable feasible region Q of biological interest 


1 TT 
ve {(,l,7, TV) ERS:L+L,< re Z “Tt 


We note that the model describes a population and therefore it is very important to prove that 
all the state variables ( L,L,;,T,T; and V) are non-negative for all time. More precisely if the system of 
equation (3) has non-negative initial data, then the solution will remain inside ( for all time t > 0, i-e., 
the set (©. is positively invariant. 











Theorem 2.2: The region ©. is positively invariant. 





Proof: In our model we have three sub populations namely; Langerhans cells, T cells and virus. 
Considering the total population of Langerhans cells only we have 





Liop =L +L; 
Hence hese _ Oy. aia 
dt dt dt 
= (a, — B,VL — uL) + (B,VL — wl, — 6 L;) 
=m, —-MWL+L;)—6L, 
Ltot 
= Ty — UL, — 6 Ly 
< 1,—-LULtt,  sinced,L; => 0 
a t+ HM Ltot ST (4) 
After multiplying both sides of equation (4) by its integrating factor e“‘, we get 
dL 
(ert) (2) + Ce" Qu Leoe S mek" 
(Leore”")' = met? (5) 


Then integrate both sides of equation (5) with respect to t. 
Eve!" = 7 e#® 4+ c, , where c, is a constant of integration. 
hope = = +¢,e-"* 

At t=0 wehave Lio¢ = Ltot (0) which implies that c, = Lyo,(0) — = 
Using this value of c, we get: 

Leot S “1+ (Lot (0) -_ “Hew 
As the number of Langerhans cells increases (i.e. t ~ 00) we have 

Ty, 
lim Ltot S cz = My 

Similarly the total population of T cells is given by T,,, = T + T; 


Thus we have 
AT +54 : dT aT, 


dt dt dt 
= (a7 — oT — B,VT — B3L,T — ByTT) + (BVT + B3l,T + ByT;T — 0 T; — pT) 
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=mt,-o0T—oT,-—pT,=t,—o(T +T,)—-pT, 
T; 
tot 


dT, : 
Ts < 17 — O Trot, Since p ,T,; => 0 


AT rot 
dt + O Trot < Tr (6) 


Multiply both sides of equation (6) by its integrating factor e? ‘. 


dT; 
(e° )—S + (€7)o Trot S (07) tr 


(Trore? )! S mre? * (7) 
Then integrate both sides of equation (7) with respect to t. 
Tie = ao +c), where Cc, is a constant of integration. 


Trot S = + c,e~% 
Att = 0 we have 7,5; = T;ot(0) so that we have cz = T;o,(0) — = Therefore we have 
wr TOUS. acs 
Trot S - ay (Teoe(0) = —) ewe 
As the number of T cells increases (i.e. t — 0 ) we have 
ry 
<— = 
lim Trot S . M2 


Now, since 0 < Leg¢ S m, and 0S Tio¢ S Mz, we have L; < m, andT; < m3. 
av 


Thus the equations of the rate of change of virus sub population a N6é L,; + MpT; — aV canbe 
rewritten as: = < Nédm,+Mpm,—-av 
a +aV < Nédm, + Mpm, (8) 


Again to solve the first order linear differential inequality (8), we need to find it's integrating factor and 

solve it as we solve for L;,, and T;o;. Therefore, its solution is given as 

wom +épm 4 ( “om + Meme 
a a 


—at 


Vit)< Oe 


: ‘ N 56m,+Mp m 
As time progress we have lim V(t) < oo s 
—>0O 


M3 

Now, since all solutions of Lio, Ti9_ and V are bounded by Q where Q = max{m,,m,,m3}, we have 
L,L,,T,T,V < Q 

Therefore, all the solutions remain in the region 0 and thus we have proved that the region 1 is 


positively invariant. 
Theorem 2.3: All the solutions of system of equation (3) are bounded. 
Proof: In proof of Theorem 2.2 we have shown that 


N 6m,+Mp m2 
ee :=™M3 


lim Liot S = = my, lim Trot S “tT = m, and lim V(t) < 
Then, we have 

Os Lem,0SL1,5m,0S5 TS m,0<S T,<S m, and0< VS m3. 
Therefore the Theorem is proved. 
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2.2.3. Steady States and Stability 
In this section, we are going to determine the equilibrium point and stability properties of the system 
of equation (3). 
Proposition 2.4: The reasonable initial conditions of the system (3) are: 


1 TC 
Lo = L(0) = sp bto = 1L,(0) = 0,7) = T(0) = mals = T,(0) = 0,V% =V(0) =0 











Proof: In the absence of virus, since L; = 0,7, =0 andV = 0, the rates of change of population 
densities of the system of equation (3) can be reduced to 


= = 1, — UL 

aT (9) 
— =, —oT 

dt 


Now, to get the steady state value of equation (3) in the absence of virus, we equate the right hand side 


of equation (9) to 0 as follows. 
m1, — wL = 0 and 2; — oT = Owhich implies L = = and T = = 


Thus, the reasonable initial conditions for infection by free virus of equation (3) can be 


Ty Ty 
L(0) = Ly = qb (0) = 0,T(0) == —,T,(0) = 0,V(0) =% =0 
and the equilibrium point of uninfected system is & 0, i 0,0). 
Endemic Equilibrium (EE) 


In the presence of virus we can find the critical point, endemic equilibrium point, of the system of 
equation (3) as follows: 


Suppose that 

f,(L,L,,7,T,,V) =m, — ByVL — ph (10) 
f2(L, L,,T,T),V) = ByVL — pL, — 6L, (11) 
f3(L,L,,T,T,,V) = tp — oT — B,VT — B3L,T — ByT)T (12) 
fa(L,L,,T,T;,V) = BVT + B3L;T + ByT;T — oT, — pT, (13) 

fs(L, L,,T,T;,V) = N6L; + MpT; — aV (14) 


Then we solve for f, = 0, f = 0, f; = 0, fg = 0 and f; = 0 simultaneously. 
From equation (10) and equation (11) we have 
m —(ut+o)L, 
BV + yu 7 Biv 





which can be rearranged as 
1, PyV — (u+ d)L; BV = (ut 6)pL, (15) 
Again, from equation (12) and equation (13), we have 
Tr _ (o + p)T; 
o+ BV + B3L, + ByT, B2V + Bgl, + B4T, 
After rearranging it, we obtain 
TB2V — (a + p)B2TV = (o + p)T (0 + B3L, + BT) — Wp Bh3LT, — Wr B4T; (16) 
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From equation (14, we have 
v= wenn (17) 
Substituting equation (17) in equation (15) and bringing like term together gives: 
Nr,B,6 Mt NB,6 (ut+6 M +6 
i LP ie 2 | By . 1 +L,T, | Pets | 


a 





— nu +8)] +7, 


ame = (wu + 6),b = eL c= “ees and d = ee 
cL? + dL,T, (18) 

Similarly, insert equation (17) into (16) and collect like term together. We obtain 
eL, + fT, = gL + hT?. (19) 


Where ¢ =~ 4 nf, f =~? — ola +p) + mrhy, g =X 4 B(5+p) and h= 


MB2p (6+p) 
a + (6+ p)By. 


Now, after rearranging and combining together the equations (18) and (19) we will obtain equation (20) 


If we leta = we will have, aL, + bT; = 


below. 
T/[gch* — dhg*] + T?|2hg(de — cf) — hg(de — ga) — g*(gb — df)] + T?[gcf? + 
(de — ga)(he + fg) — dhe? + 2eg(gb — df)| + T,[fe(de — ga) —e*(gb—df)|=0 (20) 


Suppose 
k = gch? —dhg?, | = 2hg(de — cf) — hg(de — ga) — g*(gb —df), m= gcf? + (de — ga)(he + 
fg) — dhe? + 2eg(gb — df) andn = fe(de — ga) — e*(gb — df). 


Then, equation (20) become 
kT? + IT? + mT? + nT, = 0 
but T; = 0 is the trivial solution. Therefore, we will have the cubic polynomial 
kT? +177 + mT, +n=0. 
We can solve this cubic polynomial using mat lab and only take the real and positive T; and neglect the 
rest. So, we get 





1/3 
(-27k8n+9k7 lm-2k613 +3327 k16n2—18k19 lmn+4k m3 +4k1413 n—k14 12m?) 


Le 3x 21/33 7 


1 
= Sp 472 
23(2187k5c—729k*1?) (21) 





1 
2187k3(-27k®n+9k7 lm-2kolm+3¥3V 27k n2—18k 5 mnt 4k5m3+4k4 nk 1? m? ) 
a l 
3k 
Using this solutions of T;, we can solve for L,L,,T and V and obtain 
T(hT]— NS Lj+Mp T 1 o+p)T 
| res PD y= i eat = L = (o+p)T] (22) 


e-gT| a BVH" BoV+B3Lr+BaTy 
Generally, we have seen that the system of equation (3) has two steady states: 


The uninfected steady states (VFE) (Ey = *,0,"2,0,0) and the infected steady state (EE) E* = 





(L*, L},T*,T7,V*) where L*,L,T*,T/,V* are values of L,L,,T,T;,V given in equation (21) and (22) 
above. 
2.2.4. Local Stability of Virus Free Equilibrium (VFE) 
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To discuss the local stability of virus free equilibrium, (Ey = a 0, = 0,0), we consider the linearized 


system of equation (3) at EF). 

Now we find the Jacobean matrix from equation (10)-(14) as follows. 

df, 2h af, %h af, 
aL OL; aT OT; av 
Of, Of2 Of, Of2 Of, 
aL aL, aT aT, a 
Of Of: af, oh af, 
OL OL, or OT; OV 
Of, Of Of, Os Of 
aL OL; OT OT, av 
Of; Ofs 9fs Ofs Ofs 
@L aL, OT OT, OV 


Jy Li, T, Ti, V) = 



































—BV-p 0 0 0 —BiL 
BiV —(ut+d) 0 0 pyL 
= 0 —f3T —0 — BV — B3L, — BT; —$,T —B2T (23) 
0 B3T B2V + B3L, + BgT; BT —(O +p) BT 
0 NO 0 Mp —a 
The Jacobean matrix of virus free equilibrium Ep is 
—B 1 
0 0 ByT, 
Hy 0 Lu 
0 —-(ut+6 1 
: 7 ) 0 0 By, 
Ty . Ur 0 — ‘a 
J (, 0,—, 0,0) 7 oO sy —Bytr — fat r 
o oO 
1 18 
0 B31 0 Baty S549) Bot r 
oO oO Oo 
0 NO 0 Mp —a 
Then, the linearized system is: 
du —fity 
dt — 0 0 0 7 
dv u 4 & 
“ vil oe | 0 fe 
Wa ee wi —| 0 re - - mw 
dt =1(% ae ,0,0) | e —o = Borer e 
dx oO oO 
dt y Barr Batr B2ur y 
dy ° o ° o (7 +p) o 
dt 0 Né 0 Mp —a 
=A 
If we let 
_ = _ Bi 1, 7 _ 7 B3Mr 7 - B3Mr 
Q41 =U, Qy5 = Ags = ——, A227 =U+ 6 , Azz = Aqg = ——, 033 = 0,034 = , 
L oO o 
BoM r Baty 
Q35 = 4g =~ A4g = 9 FP — ,As2 = NO ,ds4 = Mp ,as5 = a 
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, we will have 


A141 0 0 0 —Ays5 
0 —A22 0 0 ais 
A= 0 432 433, A340 435 
0 a32 0 Ag, Q35 
0 as2 0 As4  —As5 
Then, the Eigen polynomial of the coefficient matrix A is obtained as follows. 
—a,,—-A 0 0 0 45 
0 —Az2 —A 0 0 ays 
|A-—AlI=0 © 0 —A32 —A33 —A —A34 —Ag5 =) 
0 a32 0 —a44—A a35 
0 as2 0 as4 =as5 A 
A° + b,A* + b,A3 + b3d? + byA + bs = 0 (24) 


Where 
by = Qyy + yz + G33 + Ogg + Ass 
by = A441 (22 + Ag3 + Aggy + Ass) + 22033 + (Azz + G33) (A4q + Ass) + Aqqs5 + As4Q35 — Ay5A52 
bz = Q44422433 + Ay1(Az2 + A33)(Aqq + Ass) + Aq1Aq4Ass + G41 A54A35 — Aq4Qy5A52 + Az2A33(Agg 
+ dss ) + 44 Ass (Azz + 433) + AsqQ35(A22 + A33) — 45432454 + Ay5A52(A44 — 233) 
D4 = 041422033 (Aq4 + Ass) + A41A44A55 (G22 + A33) + 041454035 (Az2 + G33) — A41A45A32A54 
+ 041 A45A52 (G44 — A33) + Az2033A44As5 + A2203354035 — Ay5A32033As54 + Ay5A572033A44 
bs = Q41422033(A44As5 + Asqd35) + Qq1415433(A52A44 — A32A54) 


By Routh-Hurwitz criterion, it follows that all eigenvalues of equation (24) have negative real parts if 
and only if it satisfies: 


by SO torrets. 0 (25) 


by byb2—b3 
Note that equation (25) was obtained in section (1.7). If condition (25) holds, the virus free equilibrium 


Ey = (= O 0,0) is locally asymptotically stable. 


2.2.5. Stability of Endemic Equilibrium 
In similar fashion to the local stability of virus free equilibrium, the local stability of endemic 
equilibrium, E*(L*, L;,T*,T;, V*), is obtained by linearizing system of equation (3) at E*. Using equation 
(23), the Jacobian matrix of E* is 


_pvi-p 0 0 0 BL" 
BV" —(u + 6) 0 0 B,L* 
J(E*) = 0 —B3T*  —o0— B2V* — B3L, — ByT; —B4T* —BT* 
0 B3T* B2V* + Bgl, + B4Ty B,T°—(o+p) B,T* 
0 N6 0 Mp —a 


The linearized system of equation (3) is 
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dur 
dt 
av! —p,V* = Ul O 0 0 —B,L" u! 
ie Biv" Ce) 0 0 B,L* v! 
ae 0 —B3T* —o — B2V* — B3L; — ByTy —B,T* —B,T* || w’ 
axr 0 ite B.V* + B3l* + BsT? BxT*-—(o+p) 8,T" x 
dt 0 Né 0 Mp —a y' 
will =B OO 
dt 
(26) 
Then 
Sh 0 0 0 —bi5 
bo —b22 — Xr 0 0 bis 
IB-A|=0e0 0 —b32 —b33 —A 34 —b35 | = 0. 
0 b32 bas Dyas -A Daa 
0 bs2 0 bs4 —bss ca Xr 


Where by; = B,V* +p, bys = ByL*, bai = ByV", bez =H+ 5, bes = bys, b32 = B3T*,b33 = 0 + B2V* + 
BsL, + ByTy, b3z = B4T*, b35 = B2T*, byz = b32,by3 = Bz V+ B3Ly + B4Ty, bag = B4T* — (o+ 
P), bas = gs, bs2 = N6,bs4 = Mp, bss = a 
This implies 

P+ Gd? +A yd? + GA oe =H 0 (27) 
Where 
Cy = b33 + bss + bg2 + by, — bag 
Cz = by2(b33 + bss — bag) + b33(bs5 — baa) — bag bss — D54b35 + b34b43 — by5bs52 + by1(b22 + b33 — 
bag + bss) C3 = b22b33 (bss — bag) + b22(b34ba5 — Dagbs5 — 05435) — 033 (Daabss + bs4b35) + 
b34b430s5 — by5b33bs2 — b32bs4by5 + by5bs52b44 — 041022 (bag — 533 — Ds5) — b114b33(ba4 — bss) — 
by bygbss — by bs4b35 + by, b34b43 — by bys bs2 + boy bisbs2 
C4 = b22b34b43bs55 — b22b33(b44bss + bs4b_35) — bys5b43(b34bs2 — b32bs4) + bi5b33(bs2ba4 — b32bs4) 

+ by1b22b33 (bss — Dag + by1b22(b34b43 — bagbss — D54b35) — b11033 (Daa bss + bs4b35) 

+ by1b34b43bss5 — by1b15b33bs2 — by1b32bs4by5 + 41 Ds2b44 + b21b15(b33bs2 + b33Ds4 

— bszb4q) 
Cs = by1b22(b34b43bs55 — b33(Daabss + bs4b35)) — by1b15b43 (634052 — 32054) — b11b15b33 (b32bs4 

— bs2b44) + by5b21(b43 (b34bs2 — b32bs4) + b33(b32bs4 — bs2b_44)) 
By Routh-Hurwitz criterion described in section (1.7), it follows that all eigenvalues of equation (27) 
have negative real parts if and only if it satisfies: 
Cc, > 0,c, > 0 a2 36 — > 0,c3 — calerca ~ €3) = $5) 

Cy C1C2 — C3 
7 C5 (CyC2 — C3)? 
C3(CyC2 — C3) — Cy (C1C4 — Cs) 

Therefore, the endemic equilibrium point E* = (L*,L’,T*,T;,V*) is locally asymptotically stable if 
condition (28) holds. 


> 0,C€4C4 — Cs (28) 


2.2.6. Global Stability 
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The global stability of obtained equilibrium can be analyzed by transforming the system of equations 
into a linear system and then choosing a suitable Lyapunov function to analyze each equilibrium point. 
By letting 

L=L*+x,,L,=L),+%2,T =T*+%3,T, = T; +x4,V =V* +xs5 where x1,X2,%3,xX, and x, are small 
perturbations about L*, L’,T*,T; and V* respectively, the system of equation (3) is turned into a linear 
system which is of the form x, = J(E))x;, where J(E) is the Jacobian matrix of equations (3) and i = 
1,...,5. Thus, the linear system of equation (3) is 


dx, . : 

“dt = —(B,V* + wx, — B,L*xs 

dx 

4 = BV" x, — (u + 6)x2 + ee 

dx 

a = —B3T*x, —(o + B2V* + Bgl, + ByT/)x3 — BT *X4 — B2T*Xs5 

“ak = B3T*x2 + (B2.V* + Bgl, + ByT/)x3 + (BT =e cp p))x4 + pT *xs 
dx. 

aH NOx +Mpx,-Q@Xx-s 





Theorem: The equilibrium E (L*, 0,T*, 0,0) = & 0, = 0,0) is globally stable provided that 
O(NL* — 1) < wand p(MT* - 1) So. 





Proof: 

We define a Lyapunov function as Z(%1,%2,%3,X4,X5) = ae ae at ae xs, Where L*, T* are 
components of the equilibrium point E,(L*,0,T*, 0,0) = Ge 0, s 0,0). Since 

(x4, X2,%3,X4,X5) © R8,Z (x4, Xz, X3,X4,Xs5) is a positive definite function. Differentiating Z with respect 


to time t we get 
dZ OdZdx, OZdx, OdZdx, OdZdx, OdZ dx, 


dt 0x, dt | Ox, dt | 0x3 dt |x, dt Ax; dt 


= [-(B,V* + w)x, — ByL*xs] . + [B,V*x1 — (u+ 5) x2 + ByL*x*] 7 


i if 
1 
+ [-B3T*x2 —(o + B2V* + Bgl) + ByT/)x3 — ByT*X4 — B2T *xs] T 


+ [BsT*x2 + (B,V* + BL} + ByT7)x3 + (BiT* — (0 + p)) x4 + BoT*xs] - 
+ [N6 x. +Mp x, —-a xs] 

= ae = (“- No) x2 ae = (<P - Mp) x, —aXxs 
Choosing 








— Mp, (29) 
We have 





dZ ul U+o o o+p 
=f (CE— te) for - (GRE Ma)na ma <0 


Then, equation (29) can be rearranged as 6 (NL* —1) <u, p(MT* —1) < oso Ey = (L*,0,T*, 0,0) _ is 
Lyapunov stable. Therefore, the theorem is proved. 
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3. NUMERICAL SIMULATION 
In previous sections we have investigated analytically the properties of proposed model. Now, we shall 
investigate numerically the results of proposed model by means of numerical simulation. We solve our 
system of equation (3) by using an iterative method with Runge-Kutta fourth order scheme. 


3.1. Implementation of Runge-Kutta Fourth Order Method for Numerical Solution 
Firstly, we give a brief description of the Runge-Kutta method of order four (RK4) for the system of 
equations (3). We first develop a sets of k's used to make predictions of the dependent variable at the 
midpoint of the interval. These are then employed to make predictions at the end of the interval that 
are used to develop the values of k's at the end interval k,. Finally, the k's values are combined into a 
set of increment functions and brought back to the beginning to make the final prediction [1]. The 
following illustrates these approaches. 

We start with initial conditions L(to) = Lo, L; (to) = Lry,T (to) = To, T1(to) = Tig, V (to) = Vo- 
We assume that the values L;,L,,, Tj,T;, and V; has been computed. Now we calculate 
Lisa Li. Tita T1,,, and Viz by Runge-Kutta fourth order method as: 


Ligg = Lit (ka + 2(K21 + k31) + kar) 
Lies = Li, a = (ka + 2(k22 + k32) + Raz) 
Tis. = T; +2 (kis + 2(k23 + k33) + kas) 
ce = qT, + (ka + 2(ko4 + k34) + Raq) 


1 
Vinw = Vit G (kis + 2(kos5 + kgs) + Kas) 


With, 
kyj = hfj (to, Lo, Lig To, Ti,» Vo) 

kaj = Afi(to +=, L; +H, Ly, +27, +20, +H +29 

kaj =Afi(to +=,L, +, L, +27, +27, +H Yt) 
kaj = hfj(to th, Li + k31, Ly, + k35,T; + k33,Ty, + k34,V; + k35) where f;'s are equations given in (10- 
14) for j = 1,...,5. 
This gives us the next approximate values of L,L,;,T, T; and V. Then t is set to tj + A and the values of 
L,L,,T,T, and V are iterated with the above formula. 

3.2. Values of Parameters and Initial Variables 

Choosing variable and parameter values characteristic at the cellular level is difficult. In our model, if 
measurements have been attempted, the values taken may not be as accurate as we need for 
quantitative predictions. Thus one role of modeling is to point out where further quantitative 
measurements can improve our understanding of the AIDS disease process. For example, the number 
of CD4+ T cells in the peripheral blood is approximately 1000/mm?, although it fluctuates from time to 
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time depending on the total lymphocyte count [13]. As it is common in the clinical literature, we shall 
report all cell numbers per cubic milliliter. For our simulations, we declare the numerical parameters 
and variable by using similar procedure as for T cells above. We list these values in Table (2) and using 


these values we get Figure (2). 


Parameters Approximated Values 

Initial population size of healthy Langerhans cells (L) 800mm73 

Initial population size of infected Langerhans cells (L,) 0.033mm-3 

Initial healthy T-cells (T) population size 1000mm-3 

Initial population size of infected T-cells (T;) 0.045mm-3 

Initial HIV population size (V) 5mm~3 
The constant rate of production of susceptible Langerhans cells 19mm“? 
from the bone marrow (7;) 

Infection rate of susceptible Langerhans cells by free virus 0.002 day"? 

(B:}) 

Natural death rate of Langerhans cells () 0.002 day~+ 
Death rate of L, due to virus (6 ) 0.007 day~+ 

The constant rate of production of susceptible T-cells from the 15mm? 

bone marrow (27) 

Natural death rate of T cells (a) 0.001 day~+ 

Infection rate of T-cells by free virus (£2) 0.0002 day? 

Infection rate of T-cells by infected Langerhans cells (B3) 0.0004 day? 

Infection rate of T-cells by infected T-cells (f,) 0.003 day~+ 

Death rate of infected T-cells due to virus (¢ ) 0.008 day~+ 

Number of the virus particles assumed to produced by the 2aay > 

infected Langerhans cells during its life time including any of 

its daughter cells(N) 

Number of the virus particles assumed to produced by the 4 day~+ 

infected T-cells during its life time including any of its 

daughter cells(M) 

virus death rate (a) 0.8 day? 


Table 2: The approximated values of variable and parameters declared in the numerical simulation 


Figure 2 illustrates the dynamics of the disease for Healthy Langerhans cells, infected Langerhans cells, 
Healthy T-cells, infected T-cells and Free Virus with increasing time. The population of susceptible 
Langerhans cells decreases as time increases due to the presence of virus as shown in Figure 2(a). As 
we can see from Figure 2(b), (d) and (e), the Population of infected Langerhans cells, Infected T-cells 
and Free Viruses increases when time increases. On Figure 2(c), we observe that the population size of 
susceptible T-cells increases starting from its initial population size since the healthy T-cells defense 
viral infection, but as time increases it start to decrease. This is due to structure HIV virus. That means, 
HIV is an enveloped retrovirus. When it leaves a host cell it takes a part of that cell and duplicates itself 
along with host cells. Therefore, it starts to multiply itself within a short period of time including its 
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daughter cells so that it wins healthy T-cells. For the sake of easy comparison, we can observe all of this 
dynamical interaction on one plane in Figure 2(f). 


a) (0) (c) 
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Figure 2: Dynamics of the model showing the virus progression in an increasing time t 

4. CONCLUSION 
Infectious diseases are a threat for our modern society due to the increased population density on Earth, 
so being able to understand and predict the dynamics of infectious diseases is extremely important. 
In this research work, we have developed a model used to explain the interaction between host cells 
and HIV infection. The model was developed in accordance with previous models in the literature 
review, and in accordance of the biologically reasonable assumptions and parameters. It was shown 
that, the model has a biologically meaningful region in which we had carried out all our model analysis. 
Also, it was shown that the solution of the model systems exists and unique, positively invariant and 
non negative. The stability properties of both virus free and endemic equilibrium were determined; as 
the stability nature of solutions will assist us to determine the extent to which the disease will disappear 
from the population. In order to insights into the HIV/AIDS dynamics at cellular level, the numerical 
simulations of the model are carried out. 
In addition, in this work we have improved the understanding of the interaction between Langerhans 
cells and HIV infection by explaining the most crucial point take place in developing our model system. 
Furthermore, by using the language of mathematics we have provided the tangible understanding on 
the viral pathogenesis and its life-span for the immunologist, for the researchers and for the clinicians 
to propose the new powerful tools on the stimulation of the immune system in order to increase its 
efficiency in the struggle against antigen invasion. 


Generally, this study assists us to consolidate our knowledge on the dynamics of HIV - host cells 
interaction and it concretely simplifies the complexity to understand the interplay between the virus 
and host cells. At the future, it will be helpful to obtain good estimates for introduced parameter values 
especially, the reactivation rate of individually infected cells. Our work may require additional 
information to improve these findings by using future research studies, extensions, modifications and 
analysis of the model. Moreover, based on the parameter values, the equilibrium point may be created 
or destroyed or their stability may change. To understand these changes the model requires bifurcation 
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analysis by varying the values of these parameters. Thus, as a future prospect, it would be important to 
distinguish more clearly, on these analyses and provide a more reasonable fact on the interplay between 
Langerhans cells and HIV infections. 
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